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Abstract 

We analyze the structure of two dimensional disordered cellular systems generated by exten- 
sive computer simulations. These cellular structures are studied as topological trees rooted on 
a central cell or as closed shells arranged concentrically around a germ cell. We single out the 
most significant parameters that characterize statistically the organization of these patterns. 
Universality and specificity in disordered cellular structures are discussed. 

PACS numbers: 82.70Rr, 02.50.-r, 05.07.Ln 

1 Introduction 

Cellular structures, space-filling disordered patterns, are widespread in nature [[l], ||, In an ordered 
structure, one must constrain the elementary bricks (or group of them) to satisfy the local rotational 
symmetry compatible with translational invariance. By contrast, a disordered structure is free 
of any local symmetry constraint, and only subjected to the inescapable, topological condition 
of partitioning space. In two dimensions, these structures (froths) are partitions of the plane by 
irregular polygons. Disorder, or absence of symmetry, imposes minimal incidence numbers (3 edges 
incident on a vertex, 2 faces incident on a edge). The Four-Corner boundary between the States 
of Arizona, Utah, Colorado and New Mexico is a decision made in Washington; it has nothing 
to do with population dynamics or with agricultural efficiency. It is neither topologically, nor 
structurally stable: a small deformation splits it into two topologically stable three-State vertices, 
found everywhere else on Earth. 

The set of all the possible configurations of a froth can be explored by elementary, local topological 
transformations Q which, in two dimensions, are: Tl, exchange of neighbours (for example, a flip 
between the stable alternatives of UT-NM or of AZ-CO as neighbours, about the unstable four-corner 
boundary); T2, disappearance (or apparition) of a triangular cell into a single vertex. Two alterna- 
tive, equivalent transformations are the cell-division (called mitosis in biology) and cell-coalescence. 
It can be easily seen (fig.[l]) that local, specific combinations of division and coalescence generate 
the two elementary transformation Tl and T2, and vice versa. These elementary, local topological 
transformations are observed in all natural foams, soap froths, metallurgical grain aggregates, but 
also in biological tissues (epidermis) H or in the cellular networks produced by Benard-Marangoni 
convection, whether in the laboratory or on the surface of the sun (solar granulation) || . 

In infinite two dimensional (2D) Euclidean froths (or froths with periodic boundary conditions), 
the average number of neighbours per cell is fixed by the Euler relation to be equal to 6 |2j . Given a 
cell with n edges, the quantity q = 6 — n is called its topological charge. The sum of the charges of 
the cells over the whole infinite Euclidean froth must be equal to zero. The topological charge cannot 
be generated or destroyed. It is just shuffled in between cells by the topological transformations. 

In this paper, extensive computer simulation are performed to obtain disordered systems with 
more than 10 5 cells generated from an hexagonal lattice and applying at random either only Tl, or 
Tl and T2, or cell-division and coalescence transformations. In these simulated froths, we forbid 
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Figure 1: A froth is a cellular network generated or transformed by successive applications of two 
local elementary transformations: Tl (neighbours exchange), and T2 (triangular) cell disappearance 
(a). These two transformations can also be regarded as combinations (c) of cell division and cell 
coalescence (b). 

configurations with cells which are neighbours of themselves and cells with two edges. The typical 
resulting structures are shown in fig.^j. 

The structural organization of these computer-generated froths are studied in terms of the relative 
distances between cells. Two different topological distances can be introduced: the bond distance j6| 

and the shell distance j^, [)| [l0|, |T1J (section^). These two distances are different: the bond distance 
is defined on the graph of the froth; the shell distance is best defined on its dual triangulation. They 
are associated with the number of edges (or vertices) of the shortest paths between two cells. The 
relation between the two distances is very complex and strongly dependent on the froth organization, 
because they involve minimal paths. 

The characterization of a disordered structure is a difficult task. The absence of translation^ 
order or repetitive local configuration makes in principle necessary to know the position, size and 
shape of every cell in order to completely characterize the structure. But, disorder makes most of this 
detailed local information unimportant. The physical properties of the froth must be describable 
in terms of macroscopic statistical information. A simple and powerful way to study disordered 
cellular systems, is to analyze the structure as organized around a given arbitrary central cell. The 
information about the structure can be obtained in terms of the properties of the set of cells at 
the a given distance (bound or shell) from the central cell. We show that the classification by 
bond distances reduces the froth to a set of trees rooted on the vertices of an arbitrary central 
cell (see fig.|J). On the other hand, in terms of shell distance, the froth can be seen has a set of 
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(a) C-D 



(b) Tl 



(c) T1-T2 



Figure 2: Typical topological froths generated by: (a) Cell division and coalescence; (b) Tl trans- 
formations only; (c) Tl and T2 transformations. 

closed layers of cells arranged concentrically around the central cell (see figj^)- This classification 
in terms of distances naturally introduces a structure and a hierarchy in the froth. Moreover, a 
physical meaning is directly associated to the sets of equally-distant cells: any perturbation, signal 
or information propagates from a given cell to the whole froth through these structures and cells at 
the same distance are those reached -on average- at the same time and with equal intensity by the 
signal. They are the successive wave fronts of the signal in the froth. 



2 Topological structure of a froth 
Topological distances 

Consider two cells in a froth. These cells are connected by many paths along the edges of the 
direct graph or of its dual triangulation. Two different topological distances can be associated with 
minimal paths: 

The bond distance, is the minimum number of edges necessary to connect two cells in the direct 
graph of the frothj^. (The concepts of bond distance and of bond neighbours was first introduced 
by one of the authors who called them Tl distance and Tl neighbours, respectively 

The shell distance (or simply, topological distance (ll|), is the minimum number of edges 
connecting the centers of two cells in the dual triangulatiomg 

All the cells in the froth can be classified in terms of distances (bond & shell) with respect to a given 
"central" cell. 

We denote by F t the number of cells which are at a bond distance t and by Kj the number of cells 
which are at a shell distance j from a given central cell. These quantities are in general dependent on 
the number of edges n of the central cell. For any froth on the Euclidean plane, F t and Kj increase 
with the respective distances. But the rates of growth are characteristics of a specific cellular system 
and are good, sensitive quantities to characterize its disorder (see section |] and 

lr rhis path is a set adjacent edges ('bonds' ) connecting the two cells. 
2 The set of equidistant cells makes a closed layer around the central cell, and_successive layers form concentric 
rings. The edges separating two successive layers form a closed loop: the 'shell' [M]. 
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Figure 3: (a) The forest of n bond trees, (b) In a tree, any vertex at distance t is connected with 
one bond to a vertex at distance t — 1. Two of the three cells incident on a vertex at distance t are 
also incident on a vertex at distance t — 1. The mean vertex connectivity (zt) tends asymptotically 
to 2, which is typical of a tree structure spanning vertices distributed uniformly (see also fig.[l^). In 
(c) vertices with connectivity 0=1 and z — 3 are labelled with white and black circles respectivelly. 



The n rooted bond trees 

Associated with the bond distance, there is a spanning forest of n trees rooted on the central cell 
and reaching every vertex of the froth in the least number of steps . These trees are the minimal 
bond paths which join the vertices of the froth with the central cell. Consider a vertex of the froth 
at a large bond distance t. It has 3 neighbouring vertices, with at least one at a distance t—1. Select 
one of them and connect it with a bond. The other neighbouring vertex at t — 1 (if there is one) is 
left unconnected. Continue down in distances until a vertex of the central cell is reached. This is 
the root of the tree, which, at this stage, is simply a branch (chain) of t bonds. Then, begin again 
with another vertex at distance t, down in distances. This second i-branch either meets the first at 
some point, or it reaches another root (vertex of the central cell). Proceed successively with every 
vertex at distance t, down in distances until the branch connects another branch or is rooted into 
the central cell. Consider then the vertices at distance t — 1 which have not yet been reached by the 
branches already constructed, and repeat the operations above. By repeating the process to all the 
vertices yet unspanned at distances t, t — 1, t — 2,... .,2, 1, 0, we generate a forest of n trees, where n 
is the number of vertices (and of edges) of the central cell. These trees are rooted on the central cell 
and span all the vertices of the froth. Fig|| (a) shows a forest of trees up to distance 3. Each vertex 
in the tree has connectivity 2 < z < 3, except the root and branch ends which have connectivity 
1 (zero if the tree consists of its root vertex only, as happens, for example, if the central cell has a 
3-cell as topological neighbour). Call z t the average vertex connectivity at distance t (1 < z t < 3 
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topological defects 



Figure 4: The shell structure. Any froth can be analyzed as structured in a system of closed 
concentric layers of cells at the same shell distance with respect to the central cell. The cells (hatched 
in the figure) in layer (j) with neighbours in layers (J — 1) and (J) only, are local topological defects 
included between the layers. 



for t > and < z$ < 1). It follows immediatelly that Vt+i = {zt — l)Vt and V\ = zqVq = Zqu. The 
number of vertices at distance t > 1 is therefore related to the number of edges of the central cell 
(n) 

t-i 

i=l 

The bond distance associated to a cell is the bond distance of its vertex nearest to the central 
cell. Consider a vertex at distance t. It has 3 incident cells. Two must also be incident on a vertex 
at distance t — 1 (they are separated by the bond linking vertex t to vertex t — 1, see fig||(b) ). 
Therefore, to each vertex at a given distance, there corresponds no more than one cell at the same 
distance, i.e. F t = v t Vt with v t < 1. In 2D Euclidean froths, a cell has 6 vertices on average, and 3 
cells are incident on any vertex. Therefore we expect v ~ 0.5 on average. 

A special case is the hexagonal lattice where there are 6 vertices and 6 cells at distance t = 1, 12 
vertices and 6 cells at t — 2, 12 vertices and 12 cells at t = 3, and, in general Vt — F t — 3(t + 1) 
for t odd and Vt — 3(i + 2) and F t = 6 for t even. We have therefore v t = 0.5 on average, with 
Vt = 1 when t is odd and v t = when t is even. This difference between even and odd distances is 
characteristic of the hexagonal lattice. When the froth is disordered, this difference between even 
and odd distances disappears and, asymptotically, the quantities Vt and F t grow linearly with the 
distance t with u t — 0.5. The growth law depends on the specific system and it is studied in section 
| and H 



The shell structure 

The structure associated with the shell distance for a typical froth is shown in fig.§. The shell 
distance foliates the froth into successive layers of cells surrounding the central cell. Most cells at a 




Figure 5: The shortest bond path crosses at least t min = 2j — 3 vertices to connect the central cell 
with a cell in the layer j. 



(shell) distance j from the central cell have neighbours at distances j — 1, j and j + 1- They belong 
to layer (j). In addition, some cells at distance j (hatched in Fig.||) have neighbours at distances 
j — 1 and j only. These cells are local defect inclusions, intercalated between layers (j — 1) and (j). 

In this paper we characterize the shell structure by the number of cells at distance j, Kj(n), 
averaged over all n-sided central cells. In the perfect honeycomb (hexagonal froth), Kj = 6j, there 
are no defects, and the successive layers have an hexagonal shape. In a random froth, the layers are 
roughly circular, but wiggly. At large distances, Kj(n) is larger than for the honeycomb, but still 
increases linearly with the distance. The rate of increase is related to the degree of disorder and to 
the number of topological defects in the system |13| . 



The combined structure 

The quantities F t and Kj are both radial properties with the angular dependence averaged out. 
In order to extract information on the angular fluctuations, it is interesting to study the interplay 
between bond and shell structures. A bond tree develops radially outward from the central cell 
along edges which are perpendicular to the shells and it broadens sideways along edges parallel to 
the shells. The distribution in the bond tree of edges parallel to the shells gives therefore a measure 
of its spread. Call N(j, t) the number of couples of cells that are simultaneously at a shell distance j 
and at a bond distance t from each other. (Clearly, J^t t) — KjN and ■ N(j, t) = F t N, with 
N the total number of cells in the system.) 

To calculate N(j,t), let one of the cells be the central cell. The other cell is at shell distance j. 
What is the bond distance t between the two cells? It is readily seen (Fig.||) that at least 2 bonds 
are necessary to go from one layer to the next. Thus, t m i n = 2j — 3 = tji for j > 2. Therefore, is 
the radial component of the bond distance. (In an hexagonal froth, t = for j > 2.) 

In disordered froths, the bond distance t is, in general, much larger than £r. Suppose, for 
instance, that a path takes in average (v) extra bonds from one layer to the next. These extra bonds 
correspond to segments of the bond tree parallel to the shell and also to extra vertices associated 
to the topological defects between the shells. From the central cell to a cell at distance j, there 
are (v)(j — 1) additional bonds and the total bond distance is therefore t = trt. + {v)(j — 1) = 
(2 + (v))j — 3 — (v). The number v of extra bonds between layers is a random variable with 
probability distribution w(v) and mean (v). When w(v) is Poissonian, the conditional probability 
P(t\j) of finding a cell at a bond distance t, given that it lies at a shell distance j, is a shifted 
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Figure 6: Distribution of the number of edges per cell in the froth generated by applying cell-division 
and coalescence transformations (C-D simulation). 

Poissonian with average t = (v)(J — 1) + 2j — 3 and variance a 2 = (v)(j — 1) (see appendix |A|). The 
number of couples of cells which are simultaneously at a shell distance j and at a bond distance t is 
therefore 

e-WV-VlMtt - l)l t - 2 J'+ 3 
N(j, t) = NK 3 P{t\j) = NK, (f _2/ +3 )! ' (2) 

The arrangement of the cells in any layer can therefore be described in terms of a single parameter 
(v) only. 

3 Simulation of disordered cellular systems 
Coalescence and Division (C-D simulation) 

The disordered cellular system is generated staring form an hexagonal lattice of No = 144400 cells, 
with periodic boundary conditions, and performing 5iVo of Tl transformations on edges chosen at 
random. Then, a cell chosen at random, is divided in two parts (mitosis) if it has more than 6 sides 
or fused by coalescence with one of its neighbours if has less than 6 sides. If it has six sides, it is 
divided with probability p<j = 1/4 or left unchanged with probability 1 — pa — 3/4. The coalescence 
is performed between the cell and one of its neighbours chosen at random. We chose a symmetric 
mode of division: the two "daughters" of an n-sided "mother" cell have both 2 + n/2 edges for n 
even, and 2 + (n + l)/2 and 2 + (n — l)/2 edges if n is odd (which one has which is assigned at 
random) . 

After iterating the coalescence-division operations for 10-ZVo times, the resulting structure (shown 
in fig.||(a)) reaches a stationary number of N — 114475 cells. The normalized distribution of the 
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Figure 7: Second moment of the edges distribution (ji^) v.s. the number of elementary topological 
transformations, in the system generated by cell-division and coalescence transformations (C-D 
simulation). The system converges toward a steady state which is independent of the starting 
configuration. If the simulation starts from a very disordered state, the system self organizes reducing 
the disorder by random local transformations. 

number of edges per cell (p(n)) is shown in fig.^J. It has a maximum at n — 6 and decreases 
exponentially for n > 6) as p(n) = peexp[— a(n — 6)] with pe = 0.37 and a = 0.9. The second 
moment of the distribution is /12 = 1.77, and the Aboav-Weaire parameter [l4|, [l5| is a = 0.64. 

This distribution is stationary under further division-coalescence transformations. It is also in- 
dependent of the initial configuration. We have obtained statistically identical froths starting from 
the hexagonal lattice (/12 = 0) or from a very disordered froth with /12 > 10. In the first case, the 
system evolves from order to disorder during the simulation. By contrast, in the second case, the 
disorder of the system decreases throughout the simulation (see fig.@). This is an example of self 
organization induced by local random transformations only. The distribution p(n) reaches a steady 
state, regardless of the parameter pd which controls the fate of 6-sided cells. However, the number 
of cells is very sensitive to pd- For pd — 1/4, once the steady state is reached, the number of cells 
in the froth remains also stable. But, when pd — 0, the number of cells in the system decreases, 
whereas, when pd — 1/2, it increases. 

Random Tl (Tl simulation) 

An hexagonal lattice of N = 100172 cells and periodic boundary conditions is disordered by flipping 
(Tl transformation) 107V edges chosen at random. (Note that the Tl transformation conserves 
the total number of cells in the system.) If the chosen edge is bounding a triangular cell, the Tl 
transformation is refused (lest it would generate a cell with two sides). The probability for a Tl to 
be refused increases during the disordering process, it reaches the value of 0.21 at the end of the 
simulation. 

A typical cellular system obtained with this simulation is shown in fig.(||b). Fig.|| shows the 
probability distribution p(n) for this froth. The second moment is fj,2 = 11.9 and the Aboav- 
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Figure 8: Distribution of the number of edges per cell in the froth generated by Tl transformations 
(Tl simulation). 

Weaire parameter is a = —0.79. The average number of sides per cell is equal to 6 (as imposed by 
the Euler theorem Q) but the distribution is an exponential with a maximum at n = 3, p{n) = 
P3 exp[— a(n — 3)] with p$ — 0.27 and a = 0.30. In this simulation, most of the cells are triangles 
because it is easier to generate a triangular cell than to make it gain an edge by random Tl. 

Random Tl and T2 (T1-T2 simulation) 

In this type of simulation, a flip (Tl transformation) is performed on any edge chosen at random, 
except if the edge is bounding a triangular cell. In this case, the triangular cell is made to disappear 
(T2 transformation). For this reason, the number of cells decreases during the simulation. We found 
that, when 5./Vo transformations are applied on an initially hexagonal lattice of Nq cells, the number 
of cells is reduced by half. But 10iVo transformations are more than enough to reduce the froth to 
3 (hexagonal) cells only. This is the minimum possible number of cells for a froth tiling a torus. 
Here each cell has 2 neighbours and 6 edges and they are arranged as shown in fig.^O). Note that, for 
this configuration, there are two families of closed paths of only 4 edges, not bounding a cell, that 
wind around the torus. No further Tl or T2 transformations can be applied on this final, minimal 
configuration without generating a cell that is neighbour of itself. 

This simulation started with an hexagonal lattice of No = 250000 cells on which we applied 5Nq 
of Tl or T2 transformations on edges chosen at random. The resulting structure had a number of 
cells N — 115149. A typical froth generated by T1-T2 transformations is shown in fig.|](c). The 
probability distribution p(n) is plotted in fig|| It is quasi-stationary (until the system is close to 
its final configuration of 6 hexagonal cell). It has a maximum in n = 4 and then decreases with 
an exponential tail p(n) — j»4exp[— a(n — 4)] with p^ = 0.84 and a = 0.53, for n > 7. The second 
moment is fj,2 = 6.33 and the Aboav- Weaire parameter is a = 0.05. 
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Figure 9: Distribution of the number of edges per cell in the froth generated by Tl and T2 trans- 
formations (T1-T2 simulation). 



4 Results and Discussion 

For each simulation, the cellular systems are characterized by the number of cells in the shell structure 
(Kj), in the bond trees (F t ) and in the combined structure (N(j,t)), at the respective distances j 
and t of the central cell. The number of vertices in the bond trees Vt and their mean coordination 
Zt are also investigated. 



Bond trees 

The averaged number of cells in the bond tree ((Ft) = ^2 n p(n)F t ) as a function of the bond distance 



t is shown in fig. 11, For the three sets of simulations, (F t ) increases linearly with the distance as 
t -> oo, with a slope 3.577 (D-C), 4.68 (Tl) and 3.112 (T1-T2). As anticipated in Section || the 
number of vertices at distance t in the bond tree is expected to be proportional to the number of 
cells at that distance, Vt — F t /v t . Indeed, we find a coefficient of proportionality that asymptotically 
tends to v t — > = 0.5. 

Asymptotically, the mean vertex connectivity z t tends to 2, which is typical of a tree structure 
spanning vertices distributed uniformly on a finite area (see fig.^(c) and fig.|l2|) . There are n such 
trees for an n-sided central cell. If the trees were not impeding each other one would expect Vt oc tn, 
asymptotically in t. But trees do impede each other. This interaction between the n trees implies 
that, after a certain distance, the leading contribution to Vt is proportional to t, regardless of n. The 
influence of the central cell is screened by the disorder jl(| , which distributes the vertices uniformly 
on the plane. This is manifest in Fig.[l3|, which shows that Vt is linear in n, with a slope independent 
off. 
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Figure 10: Two topologically equivalent views of the final configuration for the T1-T2 simulation. 
There are three cells of 6 sides tiling a torus. No further Tl or T2 transformations can be applied 
on this final, minimal configuration without generating a cell that is neighbour of itself. 
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Figure 11: Average number of cells (F t ) at distance t in the bond trees, for the three simulations 
C-D, Tl and T1-T2. The average is over all central cells. 




Figure 12: A bond tree spanning l + t(t+ l)/2, uniformly distributed vertices (here on an hexagonal 
tiling). There are t vertices at distance i, at the end of the branches, with coordination 1 (white 
circles); the t vertices on the main trunk have coordination 3 (black circles); the root has coordination 
1. The remaining t(t — 3)/2 vertices have coordination 2. Thus (z) = 2 — 0(l/t 2 ). 
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Figure 13: Number of vertices Vt in the bond trees as function of the number of edges of the central 
cell n, for several bond distances t. The data refer to the C-D simulation. Similar behaviours were 
found for the other two types of simulations. 



Shell structures 

The number of cells Kj at shell distance j (whether within layer (j) or inside defect inclusions) 
increases linearly with j after the first few layers. The dependence of Kj on the number of edges n 
of the central cell is also linear. Among all cells at distance j, the proportion belonging to defects 
is 18% for the C-D simulation, 52% and 38% for the Tl and T1-T2 simulations, stationary after 
a few layers. Fig|l4 shows the averaged (over all central cells) number of cells (Kj) at a shell 
distance j. Figure jl4 shows clearly that the rate of growth C of (Kj) is very different in the three 
types of simulations. The number of cells at a given shell distance is larger in the more disordered 
systems. Specifically, the rate is C = 19.35 for the C-D simulation (/12 = 1-77), C = 56.6 for 
the Tl simulation (fi2 = 11-9) and C = 25.7 for the T1-T2 simulation (/i2 = 6.33). These values 
are much larger than C = 6 for the hexagonal tiling (^2 — 0), or than C = 2ir expected if layers 
were circular annuli. Indeed, in disordered froths the concentric layers wiggle around the averaged 
circular annulus. This behaviour has already been observed in soap and Voronoi froths |l2| . It can 
be interpreted as an additional, negative curvature caused by the disorder, and compensated by the 
positive curvature of the defects intercalated between the layers to produce a tiling which is globally 
Euclidean. A relation for the slope C in froths free of topological defects (SSI) ans with shortest 
ranged correlations was obtained in jl2|]: C = 6 + (3 — a)/i2, with a the Aboav-Weaire parameter. 
This relation give C = 10.2, C = 51.3 and C = 24.7 for the C-D, Tl and T1-T2 simulations, 
respectivelly. The moderate agreement indicates that defects are relevant. 

The topological charge (Qj) of a cluster bounded by layer (j) (included) has been measured in 
the three types of simulation. Figure ([l5]) shows the averaged charge (Qj) as a function of the shell 
distance j (the average is over all the central cells). The topological charge (Qj) is negative and 
decreases linearly with j (similar results were found for soap and Voronoi froths (l2)). The mean 
topological charge per defect ((q de ^)) is independent of j at large j. We found (q de f) = 1.1 in the 
C-D simulation, (q def ) = 2.1 and (q def ) = 1.8 in the Tl and T1-T2 simulations, respectively. 

The negative topological charge of the cluster (Qj) is due to the wiggling of the boundary. The 
linearity shows that the amplitude of the wiggling remains constant with increasing shell distance. 
The boundary does not roughen as j increases. The negative charge of the cluster is balanced by 
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Figure 14: Average number of cells (Kj) at distance j in the shell layers, for the three simulations 
C-D, Tl and T1-T2 and for the hexagonal lattice (+ symbol). The average is over all central cells. 

the positive charges of the defects just outside its boundary. 
Combined structure 

The distribution N(j,t) of couples of cells at shell distance j and bond distance t is represented 
in fig.[l(] as a function of t for different values of j, for the Tl simulations. Beyond j = 10, the 
distribution is symmetric, peaked at the mean value i = 3.39j — 4.2, with a variance a 2 — 1.35j+3.2. 
The other two types of simulation exhibit similar behaviours. The theoretical distribution (a shifted 
Poissonian, eq.(^), derived in Appendix |A|) is also plotted in fig.^j (full line), with only one parameter 
(v) to be fitted ((v) measures the mean increase in lateral spread of the tree, due to disorder, it is 
defined in Section^.) The agreement is excellent. We obtain (v) = 0.31 (C-D simulation), (v) = 1.39 
(Tl simulation) and (v) = 0.84 (T1-T2 simulation). As expected, branches are longer in the more 
disordered systems. 

5 Conclusions 

We have simulated and analyzed very different cellular patterns which have only the name of froth in 
common: they are planar networks with minimal vertex connectivity. Several disordered froths with 
more than 10 5 cells were generated by computer simulations. The results presented in this paper 
refer to three types of simulations of cellular systems (C-D, Tl and T1-T2) generated with different 
techniques, which exhibit different degrees of disorder. The Tl and T1-T2 simulations generate very 
disordered systems where the number of edges of a cell fluctuates widely. (The sizes of the cells have 
not been relaxed or adjusted in our topological simulations). The C-D simulation generates a more 
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Figure 15: Average topological charge (Qj) v.s. j, within a cluster bounded by a layer of cells at 
shell distance j. The average is over all central cells. The topological charge of a n-sided cell is 6 — n 



homogeneous cellular pattern, somewhat akin to soap froths or other natural cellular structures. 
When C-D simulations are applied on an initially very disordered froth, the system self-organises, 
and a more homogeneous configuration with steady statistical properties is reached. In all three 
types of simulation, a stationary distribution of cells can be obtained. Moreover, the number of 
cells in the froth also remains stationary in the C-D simulation for a special value of the probability 
Pd = 1/4 of dividing 6-sided cells. 

These disordered cellular structures were characterized in terms of two distances from an arbitrary 
central cell: the bond and the shell distance. In terms of bond distances, the froth is structured in 
a forest of n bond trees, rooted on the n vertices of the central cell. They have asymptotic vertex 
coordination z,^ = 2 and the number of vertices Vt at a given distance t increases linearly both with 
n and with t. 

The number of cells Kj at a shell distance j from a given central cell, increases linearly with j. 
This linearity seems to indicate that froths remain Euclidean structures, and that the concentric 
layers of cells at the same cell distance j from the central cell, if they wiggle, do not roughen. There 
is no indication of fractal behaviour in our simulations (see |l7|] for an -artificial- example of fractal 
froth). The rate of increase (C) of the number Kj of cells with the shell distance j varies strongly 
with the type of simulation. It is larger for the more disordered systems, C = 6 for the hexagonal 
tiling, and C ~ 19, 25 and 56 for the C-D, T1-T2, and Tl simulations, respectively. The slope of 
Kj v.s. j is therefore a very good parameter to differentiate froths with various degrees of disorder. 
Moreover, these differences should be observable in the physical properties of the froth like the 
diffusion coefficient [|l8| . 

The topological charge Qj of a cluster of cells at shell distances i < j from a given central 
cell is negative, with a linear dependence on j. This shows that layers of equidistant cells wiggle 
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Figure 16: Number of couples of cells N(j, t) at a shell distance j and bond distance t, as function of 
t for several values of j. The symbols are data taken from the Tl simulation. The full lines are the 
theoretical solution Eq.(||) with (v) — 1.39. Note the linear correlation between the most probable 
value of t and j. 



16 



without roughening at larger distances, that this negative topological charge or effective curvature 
is entirely due to the outermost layer, and that it is completely balanced by the positive charges of 
the topological defects just outside the outermost layer of the cluster. 

In a given froth, the two topological structures bond trees and shell structure are intimately related. 
The combined structure was studied in terms of the number N(J, t) of cells which are simultaneously 
at a shell distance j and at a bond distance t from a given central cell. We find that the combined 
structure is well described by only one parameter (v), the mean number of extra bonds between 
adjacent layers. 



Acknowledgements: T. Aste acknowledges partial support from the European Union (TMR con- 
tract ERBFMBICT950380) . 

A Distribution of bond distances of cells all at the same shell 
distance 

Here, we estimate the quantity N(j,t). First consider the layer of cells at a shell distance j from 
the central cell. Let P(t\j) be the conditional probability of finding, among all the cells at shell 
distance j from the central cell, a cell at a bond distance t. This probability can be written in terms 
of N(j, t) as P(t\j) — N(J, t)/[J2 t N(j, t)]. Let w(i) be the probability to require 2 + i bonds in the 
path between a layer and the next. Therefore, any cell in layer j is connected, in the bond tree, to 
a nearest cell the layer j — 1 through a (shortest) path of 2 + i bonds, with probability w(i). Thus, 

t— 2j'+3 

P(t\j)= w{i)P{t-i-2\j-l) (forj>2). (3) 

t=0 

Let us now calculate P{t\l) and P(t\2). Cells sharing an edge with the central cell (j — 1) are at 
a bond distance zero. Thus, P(t\l) — 6t.o)- Cells of the second layer (j = 2) are connected to the 
nearest cell in the first layer through a path of 1 + i bonds with the same probability w(i) (i = 
labels those cells connected by one bond directly to the central cell). Thus, 

t-i 

P(t\2)=^2>w(i)P(t-i-l\l)=w(t-l) . (4) 
»=o 

This fixes the initial condition in Eq.(||). 

The weights w(i) can be in principle any discrete probability distribution. The simplest case is 
when w(i) — 8i t Q. In this case, Eq.(||) has the solution P(t\j) — S t ,2j-3- This solution has values 
different from zero only for t odd, it corresponds to a structure with minimal bond length throughout, 
i.e. an hexagonal tiling, except for the arbitrary central cell. 

A realistic probability function for the weights is the Poisson distribution of mean (v): 



With this Ansatz, Eq.(|^) has the solution 

"(t-2j + 3)! 



w(i) = . (5) 

7 ' 



nm = - — , o , — ■ (6) 
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The number of couples of cells at shell distance j and bond distance t is N(j,t) = NKjP(j\t). 
As function of t (with j as a parameter) N(j, t) is a shifted Poisson distribution with average 
t = ( v )(j — 1) + 2j — 3 and variance a 2 = (v)(j — 1)). 
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